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Abstract 

Variables adapted to the quantum dynamics of spherically symmetric models are 
introduced, which further simplify the spherically symmetric volume operator and 
allow an explicit computation of all matrix elements of the Euclidean and Lorentzian 
Hamiltonian constraints. The construction fits completely into the general scheme 
available in loop quantum gravity for the quantization of the full theory as well 
as symmetric models. This then presents a further consistency check of the whole 
scheme in inhomogeneous situations, lending further credence to the physical results 
obtained so far mainly in homogeneous models. New applications in particular of the 
spherically symmetric model in the context of black hole physics are discussed. 

1 Introduction 

Loop quantum gravity [H El El provides a candidate for a non-perturbative, background 
independent quantization of general relativity which has already led to several results 
concerning the quantum structure of space and time. There are, however, also open issues 
mainly in the context of understanding the dynamics and the semiclassical limit. While 
dealing with these problems in full generality is complicated, one can isolate particular 
aspects by looking at reduced situations where only a select class of degrees of freedom is 
considered. This class of degrees of freedom needs to be adapted to the physical situation 
of interest, which is most commonly done by employing symmetry reduction. In the case 
of loop quantum gravity, or any diffeomorphism invariant quantum theory of connections, 
there is a general scheme to introduce symmetries at the level of quantum states and basic 
operators 
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Homogeneous models 13 IHl El in the context of loop quantum cosmology [TUl HH 
IT^ are by now well-understood both from the dynamical point of view and concerning 
semiclassical properties CH CSl HE] • They have led to tests of and new insigths into 
the full theory, and resulted in many physical applications |IIl[IHllinillIll2IlE21l2nil21 

123 EEl I2Z|- However, field theory aspects which play an essential role in the full theory 
cannot be tested by restricting oneself to homogeneous cases such that a generalization to 
inhomogeneous models is needed. The simplest inhomogeneous model is the spherically 
symmetric one since it has only one physical degree of freedom j2Hl 1211 iHO] • States and 
basic operators jHI], as well as the volume operator [3^, have been derived along the lines 
of a loop quantization based on connection variables, resulting in explicit expressions in 
particular for all volume eigenvalues. Having an explicitly known volume spectrum, which 
is not available in the full theory IHS IHil IMl EZI, was one of the main ingredients 
that resulted in direct calculations in homogeneous models, and so one may expect similar 
applications of the spherically symmetric model. In particular the Hamiltonian constraint 
contains the volume operator in commutators with holonomies j3H|; which need to be 
computed in order to know the constraint equation explicitly. However, it turned out 
that in the quantization of 32j the spherically symmetric volume operator has eigenstates 
different from the triad eigenstates on which holonomies would have a simple action. This 
implies that even with an explicitly known volume spectrum commutators of the volume 
operator with holonomies are hard to compute in general. As a consequence, coefficients 
of the constraint equation would have a complicated form. 

An additional complication for setting up the Hamiltonian constraint is presented by the 
non- vanishing spin connection. As discussed in [1], one often has to split off non- vanishing 
components of the spin connection from holonomies along homogeneous directions in order 
to ensure the correct classical limit. The main observation of the present paper will be the 
fact that one can advantageously split off the connection components already at the level 
of states, not only when constructing the Hamiltonian constraint operator. We will show 
that this can be done in a way consistent with both the full theory and the treatment in 
homogenous models, and even leads to a considerable simplification of the volume operator. 
That this is possible depends sensitively on non-trivial properties of the spin connection 
and extrinsic curvature for spherically symmetric configurations, which also hold true in 
polarized cylindrical wave models. Thus, similar constructions can be done in other models, 
where the method of [32] to quantize the volume operator and to find its spectrum would 
not work or where calculations would be more complicated 130] • 

We will first discuss the spherically symmetric classical phase space and constraints in 
connection variables in Sec.|21 After computing the spin connection we will start SecElwith 
preliminary aspects of the classical limit, which motivates the introduction of a canonical 
transformation to variables better suited to a loop quantization and its investigation. The 
loop representation is then briefiy done in Sec. |31 before Sec. El presents the Hamiltonian 
constraint operator as the main part of this paper. This allows several conclusions also on 
general aspects of loop quantum gravity and physical properties to be discussed in Sec. IHI 
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2 Classical phase space 



A loop quantization of gravitational systems is based on real Aslitekar variables ^TJ 02] 
which are given by the su(2)-connection = + 'jK^ and its momentum, the densitized 
triad Ef. In the definition of A, FJ^ are spin connection components compatible with 
the triad, ii'* are extrinsic curvature components and 7 G M"*" is the Barbero-Immirzi 
parameter |^ 1^ . In spherical symmetry, one only considers connections and triads 
which are invariant under rotations up to gauge transformations, which implies the general 
form 

A = A^{x)A3dx+{Ai{x)Ai + A2{x)A2)di9+{Ai{x)A2-A2{x)Ai) sin + A3 cos (1) 
and 

E = E%x)A3sm^— + {E\x)A, + E\x)A2) sin??^ + {E\x)A2 - E\x)Ai)— (2) 

with real functions A^, Ai, A2, E^, E^ and E"^ on the radial manifold B coordinatized by 
X (see, e.g., The su(2)-matrices A/ are constant and are identical to r/ = — fo"/ or 

a rigid rotation thereof, which can be eliminated by partially fixing the Gauss constraint. 
The functions E^, E^ and E'^ on B are canonically conjugate to A^j., Ai and A2: 

Qb = [ dx{dA^ A dE^' + 2dAi A dE^ + 2dA2 A dE^) (3) 
27<-^ Jb 

with the gravitational constant G. 

These variables are subject to constraints which are obtained by inserting the invariant 
forms into the full expressions. We have the Gauss constraint 

G[\] = [ dxXiE"" + 2AiE^ - 2A2E^) ^ (4) 
Jb 

generating U(l)-gauge transformations, the diffeomorphism constraint 

D[N,^]= I dxN^{2A[E^ + 2A2E^ - A^E"") (5) 
Jb 

and the Hamiltonian constraint 
H[N] = (2G)-^ / dxAr(|E"|((Ei)2 + (e2)2))-^/2 
Jb 

X {2E%E^A'2 - E^A[) + 2A,E''{AiE^ + A2E'') + {Aj + Aj - 1)((^')' + (S')') 
-(1 + 7') {2K,E%K^E' + K2E') + {Kl + Kl){{Ey + {EY))) 
=■ -H^[N]+P[N] (7) 

where is the first (so-called Euclidean) part depending explicitly on connection com- 
ponents and P the second part depending on extrinsic curvature components (which are 
themselves functions of and Ef). 
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In jST] variables 



A^{x) := V^i(x)2 + A2(x)2, (8) 
E^^ix) := V^i(x)2 + E2(x)2 (9) 

and a{x), P{x) defined by 

A^{x) =: Aicos/9(x) + A2sin/?(x) , (10) 
A^(x) =: Aicos(a(x) + A2sin(a(x) . (11) 

for tlie internal directions 

A^(x) := (A(x)A2-A2(x)Ai)M^(x), (12) 
A^(x) := (Ei(x)A2-^'(x)Ai)/E^(x) (13) 

were introduced. Similarly, we have A;^(a;) = — Ai sin/5(x) + A2Cos/3(x) and the analogous 
expression for A^. These variables are adapted to a loop quantization in that holonomies 
along integral curves of generators of the symmetry group are of the form expy4<^A;^. 
However, E''^ is not the momentum conjugate to A^, which instead is given by 

P'^(x) := 2E^(x) cos a{x) . (14) 

Canonical coordinates are thus the conjugate pairs A^, E^; A^, P"^; (3, with 

P^(x) := 2A^{x)E'^{x) sin a(x) = A^{x)P'^{x) tan a{x) (15) 

The momenta as basic variables will directly be quantized, resulting in flux operators with 
equidistant discrete spectra. But unlike in the full theory and homogeneous models, the 
resulting quantum representation has a volume operator which does not commute with flux 
operators. This follows since volume is determined by triad components, in particular E"^ 
which is related to P*^ in a rather complicated way involving the connection component A^. 
Thus, the volume operator has eigenstates different from flux eigenstates, which makes the 
computation of commutators with holonomies more complicated [22] ■ An alternative way 
to derive volume eigenvalues in a quantization based on the variables A^j., Ai and A2 is 
being pursued in [311 HO] . 

In this paper we will considerably simplify the formalism by applying a canonical trans- 
formation such that now E"^ plays the role of a basic momentum variable. This will, of 
course, change the configuration variables which will no longer be purely connection compo- 
nents. At first sight, this seems to render the procedure unsuitable for a loop quantization 
where basic operators make use of holonomies of the connection. After a transforma- 
tion of the canonical variables, holonomies in general will be complicated functions of the 
new variables such that the new quantum representation would not be suitable for a loop 
quantization. It turns out, however, that the special form of a spherically symmetric spin 
connection and extrinsic curvature for a given triad leads to new variables which are ideally 
suited to a loop representation even from the dynamical point of view. 
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3 The role of the spin connection 



The co-triad corresponding to a densitized triad Q is given by 

e = e^A^dx + e^A^d?? + e<^A^ sin '&d(f (16) 

with 

e^ = V^ and = sgn(E^)^= . (17) 

From this form one can compute the spin connection 

e' e' 

r = -(a + pyA^dx + -^A%di^ - ^Al sin + A3 cosM^ (18) 

and the extrinsic curvature for lapse function N and shift iV^, 

K = N-\e^ - {N^e,)')A^dx + N-\e^ - N^e'^)A%dd + N-\e^ - N^e'^)Al sin dd^ . (19) 
We define the 93-components of V and K as 

r.:=-J = -^ , K,:=N-\e^~N^e'^) (20) 
which combines to = ^,JY^'+^^f^C^ since the two internal (/^-directions are orthogonal. 
3.1 Classical limit 

The form of spin connection and extrinsic curvature is important to understand classical 
regimes in which extrinsic curvature is small. Since the spin connection in general is not 
a tensor, a similar statement about the smallness of F would not have invariant meaning 
since one can always choose local coordinates such that the components F^ are arbitrarily 
small. However, in a symmetric model this may not be true for all components because 
within the model only coordinate transformations preserving the invariant form ([Q), (j21) are 
allowed. Then, some components of the spin connection can become covariant such that 
statements about their magnitude become meaningful. In the above form, for instance, it is 
clear that the x-component of the spin connection is not of definite magnitude since locally 
one can simply gauge the angle a + /? to be constant. The last term in ()18|1 . on the other 
hand, does not depend on the fields at all and is thus always of the order one. In general, 
components of the spin connection in inhomogeneous directions, just as all components in 
the full theory, do not have invariant meaning, while components along symmetry orbits 
are gauge invariant. (This follows easily from the transformation properties of a connection 
where the inhomogeneous part g~^dg only contributes to inhomogeneous directions such 
as X when the gauge transformation g is required to preserve the symmetry and is thus 
constant along symmetry orbits.) Moreover, F^ is covariant and transforms as a scalar 
because both E'^' and are densities. 
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Covariant components of the spin connection carry information about the intrinsic 
curvature of symmetry orbits. In contrast to extrinsic curvature, this is not necessarily 
small in classical regimes as can be seen for the spin connection ()18|1 specialized to the 
Schwarzschild solution. Inserting e^p = x and Cx = {I — 2M/x)~^/^ leads to F,^ = 1 + 
0{x~^) which is not small at large radii where classical properties are to be expected. 
Thus, also the Ashtekar connection component = + 'j'^K'^ is not small in this 

regime but of the order one. Consequently, angular holonomies of the form exp A^A^ 
cannot be expanded in A^ for semiclassical considerations. It is then more complicated to 
quantize classical expressions which are polynomial in A^, a prominent example being the 
Hamiltonian constraint, because they first need to be expresses through holonomies such 
that the classical expression is reproduced in semiclassical regimes. 

One can arrive at expressions expandable in connection components if one explic- 
itly subtracts the spin connection for homogeneous directions. Instead of working with 
holonomies of A^ one would then construct operators from holonomies of, in the above 
case, A^ — 1 which would be small. This procedure can in fact be described as a general 
scheme which works in all homogenous models and gives a satisfactory evolution there (i.e. 
stable in the sense of jl^ and with the correct classical limit). 

The same procedure can be followed in the spherically symmetric model for holonomies 
along symmetry orbits. However, in the variables described so far the same complication 
as with the volume operator arises: In order to subtract the angular spin connection 
components at the operator level we need to quantize the spin connection which, since it 
is related to triad components, results in a complicated expression in terms of canonical 
variables A^p and P'^. 

3.2 Canonical triad variables 

The preceding remarks have shown that many steps in the usual constructions of a loop 
quantization become much more complicated when triad components are not among the 
basic canonical variables. On the other hand, applying a canonical transformation such 
that P'^ is replaced by E'^ may lead to more complicated configuration variables which are 
no longer related to holonomies in a simple way. We will now see that the explicit form (jl8|) 
of the spin connection and ()19p of extrinsic curvature in spherical symmetry allows one to 
perform a suitable canonical transformation and at the same time facilitate the subtraction 
of spin connection components already at the state level. The subtraction in the constraint 
operator can then be done easily. There is a difference to the treatment of homogeneous 
models in jU] where the spin connection was used explicitly only at the operator level. 
However, there the subtraction could also have been done at the state level such that the 
procedure here is consistent with jU] . Moreover, we will see that inhomogeneous directions 
are not affected by the subtraction and we are also consistent with the full theory. We thus 
have crucial tests available by models which are in between homogeneous ones and the full 
theory, as discussed in more detail below after explicit constructions will be available. 

Since the momentum of Ax is already given by a triad component E^, it will be un- 
changed by our canonical transformation and we can focus on the variables A^, P"^; (3, P^. 
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Using the definitions ()14|) and ()15|) of P*^ and in the canonical Liouville form and 
trading in E'^ for P'^ results in 

P'^dA^ + P^d/3 = 2E'^cosadA^ + P^d^ 

= E^d(2A^cosa) -2E'^A^dcosa + P^d/3 

= E'^d(2A^cosa) + P''dr/. (21) 

In the last line we now have E'^ as momentum of the configuration variable 2A^ cos a, and 
the old P"^ = P^ as momentum of the angle rj := a + f3 determining the internal triad 
direction. 

As a function of the original variables, cos a looks complicated and does not seem 
to be related to holonomies. In fact, a is a function of both A^^, and the momenta P'^ and 
P^ such that it cannot be expressed as a function of holonomies in the original variables 
alone. However, the structure of p8|) and p9|) shows that there is a simple geometrical 
meaning to the new configuration variable conjugate to i?'^. Here it is important to notice 
that the internal directions along a given angular direction of a spherically symmetric F in 
()18|1 are always perpendicular to those of E (note that and A'^ are exchanged in ()18|1 
compared to 0), while the corresponding extrinsic curvature components are parallel to 
those of E. Since A is obtained by summing F and K, we can write 

A^A^ = F^A + ^K^A 

with A := A'^ and A := A|;. This implies 

A^ cos a = A^A^ ■ A = -fK^ (22) 

where the left equality is the definition of a. Thus, the new configuration variable is simply 
proportional to the extrinsic curvature component K^^ which we can view as obtained by 
subtracting the spin connection from the Ashtekar connection. Note that it is well known 
in the full theory that extrinsic curvature components are conjugate to triad components. 
But as we have seen for Ashtekar variables, this does not imply that the components 
as defined here are conjugate (while E^, E"^ would obviously be conjugate to Ai, A^ as 
well as Ki, K2). The non-trivial fact is thus that in contrast to the angular Ashtekar 
components, the angular extrinsic curvature component as configuration variable allows 
one to use triad components as momenta. As the derivation shows, this depends crucially 
on properties of the spherically symmetric spin connection and extrinsic curvature. That 
E"^ is conjugate to follows from the fact that E and K have the same internal directions, 
while the orthogonality of internal directions in F to those of E is relevant for details of 
the canonical transformation. 

3.3 Hamiltonian constraint 

In the original spherically symmetric variables the Hamiltonian constraint has Euclidean 
part 

H^[N] = -{2G)-^ [ dxN{x)\E''\-^/mAl-l)E'^ + 2 cos a A^E^'iA^ + P')- 2 sin aE'^A'^) . 
Jb 
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With the new canonical triad variables, using 

A'^ sin a = {A^ sin a)' — A^^a' cos a = T'^ — •jK^a' 

we have 

H^'IN] = -{2G)-' [ dxN{x)\E^\-'/^ ((rj + j^K^ - + 2jK^E%A, + r]') - 2E^T'^) 
Jb 

(23) 

Moreover, the Lorentzian part is simply of the form 

P[N] = -{2G)-\1 + 7^) [ dxAr(x)|E^|-^/2 ^^j^v' + 2K^K^E'') 

Jb 

which, using = ~'^x = + v' ■> combines easily with the terms already present in 
(I2S1): 

H[N] = -{2G)-^ [ dxN{x)\E--\~'/mi-rl + Kl)E'^ + 2^-'K^E^{A, + r]') + 2E--r'^) . 
Jb 

(24) 

4 Loop representation 

With the new spherically symmetric configuration variables Ax,'~fK^,T] the construction 
of the quantum representation proceeds identically to that in just by replacing ev- 
erywhere Ai^ with 'yK^ and (3 with rj. The resulting Hilbert space is then spanned by an 
orthonormal basis in terms of spin network states 



Tg,kA^) = n^xp ( like / Ax{x)dx j J]^ exp{in^'yK^{v)) exp{ikar]{v)) (25) 

with edge labels & and vertex labels /ii, G M and fc^, G Z for graphs g in the 1- 
dimensional radial manifold B. 

By definition in (jH)), A^p is always non- negative which is sufficient because a sign change 
in both components Ai and A2 can always be compensated by a gauge rotation. The ex- 
trinsic curvature component fT^, on the other hand, is measured relatively to the internal 
direction Ag in and thus both signs are possible: K^p G M. This makes the representa- 
tion technically easier to deal with, but could be done similarly with connection components 
upon defining A^ := A^sgnK^. (See also fl^ for a discussion in a homogeneous model.) 

Similarly, we immediately obtain the action of fiux operators quantizing the momenta 
E^, E'^ and P'' (using the Planck length ip = VGH): 

E%x)T,,,,, = jil'^^^^^^^^T,,,,, (26) 
/ = l^lY,^^vTg,k,^. (27) 

[ P%,fc,^ = 274 J] A:,T,,fc,^ (28) 
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where e^(a:) are the two edges (or two parts of a single edge) meeting in x. (Compared to 
the expression for P"^ in jST] there is a factor of | missing in ()27|) since E"^ is conjugate to 
^K^pjG, while P"^ was defined to be conjugate to A^jlG; see (j2T| . Note also that in some 
previous papers = SnGh was used.) Spin networks, as expected, are thus eigenstates 
of the flux operators with the crucial difference to jSI] being that now E'^ is among the 
fluxes. 

This allows us to quantize the volume V{I) = Air Jjdx^y\E^\E'^ of a region I C B 
immediately after regularizing as in resulting in the volume operator 




where ^{x) is the distribution valued operator 

E^{x)Tg^k,t. = 74 XI ^)^^'^9,k,,. ■ 

Note that, just as in ©, also E"^ is defined to be non-negative in Q. Thus, only labels 
/i^ > would be allowed. Again, it is technically easier to allow first all values /i^ G M 
and in the end require physical states to be symmetric under /i„ h-> —fi^ corresponding to 
solving a residual gauge transformation. We thus write explicitly absolute values around 
E'^^x) and /x^. The volume operator then has eigenstates with eigenvalues 

Vfc,^ = 4717^/24 \^^^\^J\\K+(^^) + K-{^)\. (30) 

V 

In contrast to [SI], the eigenvalues follow immediately and the eigenstates are identical to 
flux eigenstates. 

The Gauss constraint and the fluxes it depends on are unmodified as compared to 121] 
(except for the — sign in front of which was written by mistake in jSI]) such that 
classically it still is given by 

G[\] = I dxA(^^' + P'') ^ (31) 

JB 

and quantized to 

G[\]Tg^k,i, = 74 A(t;)(fce+(„) - + 2ky)Tg^k,fM = . (32) 

As before, this is solved directly and imposes the condition 

= ji^^'^i'") ^e~(ti)) (33) 

on gauge invariant states. Inserting this in ()25|) we eliminate the integer valued vertex labels 
and obtain the general form of gauge invariant spherically symmetric spin networks 

^s.fc.M = n^^P ( i^^'^ •^'^"^^ + ) Yi^wi^f^vjK^iv)) ■ (34) 
They depend only on the gauge invariant configuration variables + 1]' = 'jK^ and K^. 



9 



5 Hamiltonian constraint 



There is a general procedure to quantize objects such as the Euchdean part 



H^[N] = -{8nG)-^ [ d^xN{x) 



ab 



^/\detE 



of the Hamiltonian constraint in the full theory, where the curvature components 
of the Ashtekar connection are expressed via a holonomy around a small loop starting 
at X, and the factor containing triad components is obtained from a commutator between 
holonomies and the volume operator [SB] • This procedure can also be adapted to symmetric 
models, which allows detailed tests of its viability on general grounds. While the commu- 
tator part goes through almost unchanged in symmetric models (see, however, [IE]), the 
quantization of curvature components is different. For one obtains the correct components 
only in an expansion of the holonomy for a closed loop, which implies that the exponent 
appearing in this holonomy must be small. This is easily achieved in the full theory by 
making the loops small enough in coordinate length, which results automatically in the 
continuum limit in which a regulator is removed. A similar argument can be applied to 
inhomogeneous directions in symmetric models such as the radial direction in the system 
studied here. We can then use, e.g.. 



when e = f + — f is the coordinate distance between two vertices v and f + connected by 
the edge. 

For homogeneous directions, however, this argument cannot be applied since there 
is no free edge length available which could be chosen small. Homogeneous connection 
components transform as scalars and thus are implemented via 'point holonomies' ^71 0HI 
mi?] such as expijK^{v) rather than holonomies along edges. An expansion 



is then possible only in regimes where the relevant connection component A is small. 
(Alternatively, one can employ a limit 7 ^ as in JH], but it must still be shown to 
exist in an inhomogeneous construction. This may also require conditions on the basic 
algebra used.) Since the expansions are done in order to reproduce the classical limit, 
only exponents in holonomies are allowed which become small in semiclassical regimes. In 
general, there will be other regimes in which the components are not small, which will give 
rise to perturbative quantum corrections |^ I^U} 1^ IT^ 13^ 

5.1 Classical limit 

Smallness in semiclassical regimes is not guaranteed for A^, as discussed before in Sec. 13.11 
and so one would need to subtract off the spin connection from this component along the 




expiA{v) ~ 1 + iA{v) + 0(^4^) 
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lines of jH] in order to ensure the correct classical limit. This is complicated in the spheri- 
cally symmetric model since the spin connection, containing E'^, would be a complicated 
operator in a loop quantization based on configuration variables containing A^. In the 
quantization described here, and thus f ^ are simpler, but a subtraction is not even 
necessary since now is a basic configuration variable. Since the extrinsic curvature com- 
ponent is small in semiclassical regimes, we can work directly with holonomies expi'^K^ 
which are basic and also appear in spin network states (^3]) . 

When there is non- vanishing intrinsic curvature in a model, one can proceed by splitting 
the Hamiltonian constraint into two parts which will be quantized separately The 
first one just contains extrinsic curvature components besides triad components and will 
be quantized by employing holonomies around suitable loops. In our case, this part is 
classically given by the contribution 

Hk[N] = -{2G)-' [ dxN{x)\E^r'/^ {KlE'P + 2r^K^{A, + i)E^)) (35) 

J B 

to (j21I). The second part contains the curvature of the spin connection, which in our case 
results in 

Hr[N] = -{2G)-' I dxN{x)\E''\-^l^ ((1 - rJ)E^ + 2T'^E'') . (36) 
Jb 

Both parts consist of two terms containing E"^ and E^, respectively. Via a/I-E^I oc 
and a/I-E^I oc {K^,V} these terms will contain commutators hx[h~^,V] and 
hip[h~^, V], respectively, when quantized. Here and in what follows we use holonomies 



/i,(e) = exp(/A,(a;)dxA3) (37) 

e 

h4v,6) = exp{-f6K^{v)A) (38) 

h^{v,6) = exp{^6K^{v)A) (39) 



adapted to the symmetry in order to construct the operator. The edge length and the 
parameter 6 have to be chosen for the construction; their roles will become clear later on. 

The spin connection part Hr will be completed by using the expression for and 
the triad operators. Since classically the inverse of E"^ appears, one has to use techniques 
as in in order to obtain a densely defined operator. For the extrinsic curvature 

part Hx we need to construct appropriate products of holonomies in order to model loops 
resulting in the correct curvature components. Suitable products can be read off from the 
general form 

for a quantization, where we sum over vertices v and triples i^j^k of edges and use the 
lapse function smearing the constraint. Holonomies for each edge appear in the 
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commutators, multiplied with a loop holonomy hij lying in a plane spanned by the two 
other edges of the triple. In a symmetric context, this scheme takes the form 12301313 

Hoc N{v)e^^^ti{hihjhfhfhK[hK\V]) 

v,I,J,K 

where hij is constructed from holonomies available in the symmetric model, where now 
/, J, K correspond to coordinates adapted to the symmetry. In general, each holonomy hj 
for a homogeneous direction appears with a parameter 6 (which we will use but not specify 
in what follows) with 5^ being analogous to the size of the loop in hij. This parameter 
allows one to do formal expansions of functions of holonomies, and to check if the correct 
classical limit results to leading order, but this limit will also be valid for small K^, as 
expected in semiclassical regimes, even if S is of the order one. 

In spherical symmetry, {/, J, K} is {x, i), such that we have two essentially different 
contributions for K = x and K = respectively. They lead to operators 

ii{h^{v)h^{v)h^{v)-^h^{v)-^h^[h-\V]) 

and 

ii{h^hi){v+)h-^hi){v)-^h^{v)[h^{v)~\V]) 

where the x-holonomy is computed for an edge connecting the vertex f to a new one, 
f+. Using x-holonomies for small edge length e and expanding in 5 (or using small K^p 
in a semiclassical regime) we see that both terms indeed give us the expected curvature 
components. We first expand the product of holonomies 

h^h^h^'^h;^^ = exp{'ySK^A) exp{^6K^A) exp{-^SK^A) exp{-'ySK^A) 
= {l + -f6K^A + 0{6^)){l + -f6K^A + 0{6^)) 

x(l - ^6K^A + 0(52))(1 - -f6K^A + 0{6^)) 
= l + t'^'Kl + j'6'K'^A, + 0{6') 

and 

Kh4v+)h-'h4v)-' = e^^^'''exp{^6K^{v+)A{v+))e-^^-''-'exp{—f6K^{v)A{v)) 

= (1 + / A,A, + 0(e2))(l + -f6K^{v+)A{v+) + 0{6')) 

X (1 - / A^M + 0(e2))(l - ^5K^iv)Aiv) + 0(5^)) 
= 1 + ^6K^MAM - ^5K^{v)A{v) + \i5K^{v^) / A,A{v+) 

+\^5K^{v) / A,A{v) - ^H^K^{v)K^{v+)A{v)A{v+) 

which will be multiplied by the commutators hx[h~^, V] and h^{v)[h^{v)~^ , V], respectively. 
To leading order in e and S, the first commutator is proportional to A3, while the second 
one is proportional to A{v). Thus, the relevant traces to check the classical limit are 
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and 

-2tr{h^h^{v+)h-'h^{v)-'A{v)) - -2-f6K^{v+)tT{A{v+)A{v)) 

-l6K^{v+) JA^ tr(AK)A(^;)) - -fSK^{v) JA, tr(A(^;)2) 
= 'y6K^{v+) sin At] 

+ l-f6K^{v+) JA^ cos Ar] + ^-fSK^iv) JA^ 
= ^6eK^{v){A,{v) + r]\v)) + 0{5^) 

with A?7 = r]{vj^) —r]{v). With these expressions we obtain the right coefficients of and 
E'^ in Hk to ensure the correct classical hmit. 



5.2 Operator 

The preceding discussion, together with collecting numerical factors, shows that the oper- 
ator 

H[N]= ',^,^2 E ^N{v)tr(^{h,h^h,'h-'~h^h,h-'h-' (41) 

P v,a=±l 

+ (^K,ahAva)h-^Mv)-' - h^{v)K,ah^{v,)-'h-]^ + j%A{v)^ h^[h-\ V] 

+ {h^{v)h,,,K{v„)-^h-], - K,,K{v,)h-]^K{v)-^ + 27^5 / f;A(t;)^ h^[h-/ , V]^ 

corresponds to the correct classical expression where in the continuum limit coordinate 
differentials e from /?.,,. and j T'^ become the integration measure. (The sign factor a 
denotes the orientation of the radial edge running to the right (+) or left (— ) of the 
vertex. If it is not specified, it is understood to be positive without crucial changes for the 
other orientation. Moreover, in front of the derivative of the spin connection component 
indicates that it has to be integrated between v and as a result of the discretization.) To 
evaluate the action explicitly we now use exp(74A) = cos \ A + 2A sin ^A for all holonomies. 
This gives, for instance, 

+4cos(|75X^) sin^(|7(5X^)(A^ - A^) 

and 

h.[h-\V] = 1/-cos(|/^,)ycos(|/A,)-sin(i/^,)T/sin(|/A,) 

+2A3 (cos(i / A,)V sin(i / A,) - sin(i / A,)V cos(i / A,)) 
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which, when combined with the contribution where -(9 and (p are exchanged and traced, 
yields the term 



X (cos(i / A,) V sin(i / A,) - sin(i / A,)V cos(i / A,) 



as one part of the constraint operator. 
Similarly, we obtain 



-2tT \^[Kh^{v+)h^'^h^{v) ^ - h,Q{v)h^h^{v) ^h^^) h^[h^^ ,V] 

+ {K{v)h,h-\v^)h-^ - h,K{v+)h-%\v)) h^h-.W]) 
-- Acos{\-i5K^{v))iiin{\-i5K^{v+))sm{j A,-7]{v)+7]{v+)) (42) 

X (cos(i75^^M)V^sin(i75ii'^(i;)) -sin(i75ir^(t;))t>cos(i75^^M)' 



The spin connection component = — and its integrated spatial derivative V'^ = 
(v) can be quantized using flux operators, where we need to choose a discretiza- 
tion for the derivatives and be careful with the inverse of E'^ because the flux operators 
are not invertible. The latter problem can be solved in a manner by now common in loop 
quantum gravity, expressing the classical inverse as a Poisson bracket between holonomies 
and only positive powers of flux operators jSlj. The derivatives are also straightforward to 
deal with because the total expressions are scalar. We can thus write 

E^' 1 f E^M - E-{v) E^{v.)-E-{v) \ , ^ 
^-^'^ = -2E^ = -4 [ lE^ j + ^^'^ ' 

treating the two neighboring vertices symmetrically, where now all expressions in the nu- 
merators and denominators are scalar and can be quantized. Note that in this manner 
becomes non-zero only in vertices. There are several choices involved in the construction, 
choosing a quantizable form for the inverse |H7| IKH] and a discretization of the derivative, 
but qualitative aspects to be discussed in what follows are not affected. 

5.2.1 Action 

To write down the action on triad eigenstates explicitly it is convenient to split the vertex 
contribution to the operator into three parts, Hy = Hi, + He + -f^R with 

Hr/l = ^Q^^2f cosil^SK^iv)) sm{l^6K^{v±)) sin i^f A, - r]{v) + r]{v±) j A^V (43) 

depending on K,^ in v and or v^, respectively (receiving contributions from the two 
bottom lines in (j4H) for cr = + for R and cr = — for L) and 

= ^^'3^212 {sm\^6K^)A^V + ^^6^1 - fJ)A,y + ^^6T'^A^V^ + (44) 
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depending on K^p only in v (with contributions from the top hnes in ()41|) for both cr = + 
and 0" = — ) where 

Aj/ := cos(i / A,)V sin(i / A,) - sin(i / A,)V cos{\ J A,) (45) 

and 

A^V := cosi^^5K^iv))Vsm{l^5K^iv)) - sm{l^5K^iv))V cosil^5K^iv)) . (46) 

Since the expression for V in terms of /c+ and k_ only depends on the sum + k_ 
and the operator A^^V^ turns out to be diagonal on triad eigenstates, we do not need 
to distinguish between the versions integrating to f+ and f_, respectively, if v± are new 
vertices. This is different if v± already exist as vertices of the original graph which, however, 
does not crucially change coefficients. We will thus suppress this additional possibility in 
the notation. 

Acting on a vertex v, only labels of that vertex, /x, and its two neighboring ones, /x± as 
well as the connecting edge labels, k± are changed. We therefore drop all other labels in 
the notation for states \fj,-,k-, fj,,k^, fj,^), symbolically 

k- k 

which have connection representation 

{Kip, AxlfJ,-, k_, fj,, /c_|_, 



The action of the contributions to the Hamiltonian constraint then is, first, 

Hc\ll,k) = 272(773/2^2 ('^' (^'^^ + A:_ + 1| - V|A:+ + A:, - 1|) (48) 
X (|yu_, fi + 26, k+, + yu - 26, k+, 

-2(1 + 27^52(1 _ rj(/i, kMfi., /X, /c+, /i+)) 

-47^52 sgn5/2(/i) + k-\T'^{fl, k) /i, k+, 

where for T^p and its derivative we have to insert eigenvalues as functions of graph labels. 
Depending on the quantization chosen, this may require further labels beyond /c_|_ and /c_ 
written here explicitly. Similarly, the expression \J\k^ + A;_ + 1| — \/\k^ + k- — 1\ which 
occurs if v± were not vertices of the original graph can depend on other labels if those 
vertices were already present; these possibihties will be discussed below but coefficients 



exp(i/i_7K^(f_)) exp \ \ik_ J {A^ + r]')dx j exp{ifj,-fK,p{v)) 
exp (lik+ f {A^ + ri')dx \ exp(i/i+7i^'^(i;+)) . 
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here do not change the main results of the present paper. In the term containing we 
introduced the function 

^ r 1 for /X > 5/2 

sgn5/2(/^):=T(|/^ + V2|-|/^-V2|) = < 2/x/5 for - 5/2 < < (5/2 
^ [ -1 for /X < -5/2 

as it follows from A^F and also occurs in 

Hn\fi,k) = 4^(^^3/2^2 sgn^/2M + fe-l (49) 

X /X + |5, /c+ + 2, /i+ + ^5) - + ^5, k+ + 2, /i+ - ^5) 

+ k^,fi- |5, /c+ + 2, /i+ + i5) - k_,fj,- ^6, k+ + 2, /x+ - |5) 
A:_,yU + |5, k+ - 2, /i+ + |5) + + |5, A;+ - 2,/i+ - |5) 

- A;+ - 2,/i+ + i5) + A;_,/x - \5, k+ - 2,/x+ - i5)) 

and analogously for Hj^ where k- and /x_ change instead of k^ and 

It is convenient to suppress the vertex labels and write explicitly only changes in k^ with 
coefficients given by vertex operators changing only vertex labels but potentially depending 
on neighboring edge labels. To that end we introduce Co{k) := He together with 



C'R±(fc)|/x_,/x,/x+) := ^ 4^(g^3/2^2 ^S'^V2MVl^+ + fe-l + |5,/x+ + |5) 

-|//_, + i5, //+ - |5) + |//_, /X - i5, //+ + |5) 

- !//_,//- |5,/i+ - i5)) (50) 

C'l±(^) I//-, //,//+) := ^ 4^(g'y/2^2 ^g"^/2M + ^-1 (1/^- + // + |5, //+) 

- \fi- — ^S,n+ |5, + \fi- + |5, /i. — 

- - i5,/x- i5,/x+)) (51) 

and the constraint becomes schematically 

H[N]^{k) = 5]iV(^)(C'o(fc)^(...,^_,^+,...) (52) 

V 

+Cn+{k)H- ..,k-,k+ + 2,...) + Cn-{k)ip{- ■ ■ ,k-,k+ - 2, . . .) 
+C'L+(fc)^(. . . , A;_ + 2, A;+, . . .) + C'L-(fc)^/'(. . . , A;_ - 2, . . .)) . 

5.2.2 Difference equation 

Until now, the particular choice of vertices v± only affected the form of coefficients in Co 
which will not be very important for understanding the evolution scheme. Now, we assume 
that v± have already been present in the graph before acting with the constraint operator. 
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The number of labels in a state then does not change after acting and we can write the 
constraint equation equivalently as a set of coupled difference equations. 

Upon transforming to the triad representation by expanding |?/^) = J2k fii^i^^ f^) ^ 
the constraint equation H[N]\iIj) = for all becomes equivalent to the set 

C'r+(A;_, k+ - 2)V(. ..,k.,k+-2,...)+ Cn^k., k+ + 2)V(. . . , k_, k+ + 2, . . .) 
+C'l+(A;- - 2, A;+)t^(. . . , - 2, A;+, . . .) + C'l-(A;_ + 2, A;+)V(. . . , k. + 2, k+, . . .) 
+C'o(fc-,MV(---,fc-,fc+,---) = (53) 

of coupled difference equations, one for each vertex. The adjoints are taken in the vertex 
Hilbert spaces only, i.e. 

-sgn5/2(Ai + i'5)|/U-,/U + 11+ - \5) 

+ sgn5/2(/^- i5)|/w^,/i- /i+ + \5) 

- sgn5/2(^ - \5) l/i-, /i - /i+ - (54) 

So far, we used the ordering with the quantization of triad components to the right, 
and the resulting operators are thus not symmetric. After the construction we can now 
also consider other orderings, most importantly the symmetric one. For this ordering, 
the coefficient operator C-^^{k_,kj^ — 2)"^ in (|53|l is replaced by |(C'r^(A;_, fc+ — 2)''' + 
CK-{k-, k+)). From (fSUj) and one thus obtains coefficients proportional to sgn5/2(Ai ± 

\5)\J\k+ + A;_ — 2| +sgn5/2(/^)\/l^+ + ^-1 the terms multiplying ■?/'(..., fc+ — 2, .. .), 
the sign in /i ± 5/2 depending on whether is raised or lowered in the term. It is easy 
to see that these values, unlike the coefficients sgn5/2(/^) ^1^+ + ^-1 for the non-symmetric 
ordering, are non-zero for all k±ii n ^ If /i = ^fI*^, the coefficients become zero for 

kj^ + k_ = 1. The coefficients in the non-symmetric ordering, on the other hand, are zero 
if A;+ + A;_ = irrespective of the value of fi. This has consequences for the singularity 
problem as discussed briefly below and in more detail in 



5.3 Regular izat ion issues and anomalies 

The special nature of 1-dimensional graphs relevant for spherically symmetric spin network 
states makes some issues in the context of the Hamiltonian constraint more complicated 
than in the full theory. For instance, to have a well-defined operator after regulators are 
removed it is necessary that the constraint operator (jlT|l annihilates spin network states 
based on graphs without any vertices. Otherwise, there would remain infinitely many 
contributions of the same non-zero value after the continuum limit is performed. 



5.3.1 Action on special vertices 

For this issue it is sufficient to consider the commutators of the volume operator with 
holonomies appearing on the right hand side of the constraint operator. While the volume 
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operator annihilates states without vertices, this is not obvious for the second contribution 
to the commutator where we first act with a holonomy. The exphcit form of the matrix 
elements, however, shows that the commutators appearing in the constraint always anni- 
hilate states with zero vertex label /i^, such that we obtain a well-defined operator from 
the non-symmetric ordering in the continuum limit. Note that the argumentation here 
is less trivial than in the full theory, where one was able to refer to planarity of vertices 
obtained after multiplying with an edge holonomy in the commutator. Here, one has to 
use an explicit computation of the commutator which then results in the same conclusion. 

As in the full theory, however, we do not obtain a well-defined operator if we use the 
symmetric ordering before removing the regulator. Nevertheless, one can certainly take the 
continuum limit of the non-symmetric operator, which is densely defined, and symmetrize 
it by adding its adjoint. This does result in a well-defined symmetric operator, showing 
that there is still an ordering ambiguity, but it is different from the symmetric ordering 
considered before because the adjoint removes vertices. This remark also applies to the 
full theory: regularization arguments cannot fix ordering ambiguities since a reordering is 
always possible, and may even be forced for other reasons, after the construction of oper- 
ators. To fix ordering ambiguities one needs additional arguments, such as the singularity 
statements provided below. 

A closer look at the behavior in vertices reveals that the situation can be seen as related 
to planarity at least heuristically. For the conclusion that the action of the Hamiltonian 
constraint is zero if there is no vertex relies on the presence of sgn(/i^) in some terms 
which would otherwise not vanish for fi^ = 0. (Since some of the signs are replaced 
by sgn^/2(A*'y ^ in the symmetric ordering, the continuum limit does not exist in 
this case.) The absolute values, as noted earher, appear because of the residual gauge 
transformation /i^ —fJ^v, which itself is a consequence of isotropy in the spherical orbits. 
Connection components in the and <y9-directions then merge into one local kinematical 
degree of freedom, and in that sense vertices created by a single holonomy, as they appear in 
commutators, can be considered planar. This is not automatic, however, for more general 
vertices since they can be interpreted as incorporating as well as information and are 
thus not annihilated. In polarized cylindrical symmetry, where similar states and operators 
can be written, there is no such isotropy in orbits and there are two independent vertex 
labels. Planar vertices then only arise if one of the labels vanishes, in which case they are 
automatically annihilated by commutators as they always contain at least one of the two 
vertex labels factor. 

Similarly, in the full theory one can use the planar nature of newly created vertices to 
conclude that two Hamiltonian constraint operators with different lapse functions commute 
jHH] and then discuss the issue of anomalies. Again, we do not have this argument at our 
disposal since there is no notion of planarity in the reduced, 1-dimensional manifold. If 
the interpretation of planarity given above is used, there is no simple argument since now 
different holonomies, one for the vertex dependence and one for the commutator, have 
to be considered. This may indicate that the argument for anomaly freedom in the full 
theory is too simple, as it has indeed been criticized before jUTH IHT| on heuristic grounds 
of not giving local propagating degrees of freedom, and for other reasons more recently 
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in |n2]- In this context, one- dimensional models can be helpful because they provide 
two treatable situations where local physical degrees may (polarized Einstein-Rosen) or 
may not (spherical symmetry) be present (see also |^). Nevertheless, the Hamiltonian 
constraint operators appear very similar such that their properties do not obviously show 
how propagating degrees of freedom would be realized. 

5.3.2 Quantization without creating new vertices 

To understand this issue it is helpful to consider different versions of the constraint operator, 
and different regularization procedures. The main question is how the new vertices v±, 
analogous to new edges in the full theory, are chosen. The standard way is to create a new 
vertex next to v not present in the original graph. In the limit removing the regulator, v± 
approach v which allows the expansion of edge holonomies. This is the situation where one 
needs to make sure that the operator has zero action if there is no vertex present originally 
because the classical constraint regulated by a Riemann sum has contributions everywhere. 
It is also the situation where, in the full theory, the argument for anomaly freedom of |68] 
applies. 

Alternatively, v± can be chosen to be always the next vertices already present in the 
original state. This would clearly prohibit the resulting operator from being viewed as a 
regulated expression as before, since v± cannot approach f at a fixed state. The viewpoint 
here is instead that one has an operator which reproduces the classical expression only for 
suitable semiclassical states, which still is to be proven, but not on arbitrary states (which 
in general show quantum behavior or at least corrections). The first property would be 
one part of the justification for using that operator, and one condition for it to be given is 
that expansions of holonomies and exponentials are valid in semiclassical regimes, which 
we have demonstrated before. In addition, there must be other consistency conditions 
which guarantee that the solution space to the operator is big enough for sufficiently 
many semiclassical states. Here, one usually encounters the anomaly issue: If the classical 
constraint algebra is not mimicked in a certain sense by the quantum operators, the solution 
space can be too strongly restricted. In what sense this is to be ensured will have to be 
determined by understanding the physics of models. 

In the prescription for the constraint where v± are already present as vertices originally, 
the technical issue of verifying the commutator algebra is more involved but in addition 
to the fact that the symmetric ordering is here well-defined it has the advantage that the 
constraint equations can more easily be written as a system of coupled difference equations. 
When v± are created as new vertices, on the other hand, there will be new degrees of 
freedom involved after each action of the constraint which makes it more difficult to be 
formulated as a set of equations with a given number of independent variables. When 
difference equations are available, the anomaly issue is related to the question of whether 
or not all coupled difference equations are consistent with each other, i.e. whether or not 
a consistent recurrence scheme can be formulated to solve the equations from given initial 
and boundary conditions. This turns out to be possible despite of the fact that the 
absence of anomalies is not clear yet. There are thus sufficiently many solutions for this 
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version of the operator. Still, an analysis of possible anomalies is of considerable interest 
as it can, for instance, reduce the ambiguities discussed, e.g., in [Hll, but has to rely on 
explicit computations which will not be pursued here. 

5.3.3 Diffeomorphism constraint 

In the preceding paragraphs we only discussed the Hamiltonian constraint and operator 
equations it implies, but the constraint algebra also has to be tested in combination with 
the diffeomorphism constraint. As in the full theory [^3], there is no operator for the 
infinitesimal diffeomorphism constraint, but the exponentiated version corresponding to 
finite diffeomorphisms can be quantized analogously [HI]. It simply moves vertices implying 
that their position in the background manifold B has no physical meaning. Also the 
coordinate position of vertices v±, if they are created by the Hamiltonian constraint, is 
meaningless at the diffeomorphism invariant level, but creating a new vertex between two 
given ones is still meaningful at the level of diffeomorphism classes of graphs. The version 
of the constraint operator where no new vertices are created can thus be formulated equally 
at the diffeomorphism invariant or non-invariant sectors. When new vertices are created, 
on the other hand, we need to choose their positions, or work directly on diffeomorphism 
invariant classes without any such choice being necessary. 

The latter procedure, i.e. formulating Hamiltonian constraints which create new vertices 
directly at the diffeomorphism invariant level can indeed be done despite general statements 
to the contrary sometimes encountered in the literature. Those statements point out the 
fact that the lapse function appears in the smeared Hamiltonian constraint and also in 
a quantization such as ()40|) . The factors N{v) in vertex contributions of the operator are 
not diffeomorphism invariant which renders a simultaneous imposition of the Hamiltonian 
and diffeomorphism constraints impossible. 

However, after quantization each vertex contribution is a well-defined operator and has 
to be imposed independently because the N{v) are free. The smeared constraint operator 
should be seen as a collection of all the individual vertex constraint operators in a compact 
manner. In contrast to the classical case where only the smeared constraints have well- 
defined Poisson brackets, the vertex operators have well-defined commutators as a result of 
the spatial discreteness contained in a spin network state. One can thus consistently work 
with the vertex contributions for allvEB (which are zero whenever acts on a state 
which does not have f as a vertex) together with the diffeomorphism constraint. Since 
N{v) is no longer involved, there is no problem of defining the Hamiltonian constraints on 
the diffeomorphism invariant sector. 

6 Discussion 

In this paper we have constructed the Hamiltonian constraint operator for spherically 
symmetric models within a loop quantization and computed its full action. Explicit cal- 
culations were facilitated by a choice of new variables which are a mixture of connection 
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and extrinsic curvature components. This led to several simplifications as compared to 
the pure connection variables used previously in [SH 123 12^1 HO], all related to the fact 
that the volume operator simplifies. A volume operator with explicitly known spectrum is 
also available for the pure connection variables [S2], but its eigenstates are not identical to 
flux eigenstates. This fact implies complicated expressions for commutators of the volume 
operator with holonomies which appear in the Hamiltonian constraint. In particular the 
commutator [/i<^, V^] = [cos(|A;^),K] + 2[sin(|y4<^)A^(/3), with an angular holonomy is 
cumbersome since sin ^A^ has a complicated action on volume eigenstates. Moreover, in 
this framework jS does not commute with the volume operator such that there are several 
different commutators in the constraint. 

In the new variables introduced here, on the other hand, volume eigenstates are identical 
to flux eigenstates and the relevant angle rj does commute with the volume. Thus, there 
are less commutators with different action, and each of them is easy to compute. The 
relevant calculations of matrix elements of the constraint operator are no more involved 
than in homogeneous cases jHIEl; only the constraint equation itself is more complicated to 
solve or discuss since the system involves infinitely many kinematical degrees of freedom. 

6.1 Testing the full theory 

Symmetric models are introduced to test a possible full theory in simpler situations and 
to derive physical applications in a more direct way. For reliable results it is essential that 
a model is as close to a full formulation as possible in order to ensure that there are no 
artefacts from using simplifications of the model. ^ In this respect one may worry that the 
introduction of new variables here spoils the relation to the full theory since it leads to 
key simplifications in an otherwise barely treatable system. Yet, the quantization in new 
variables is in many respects closer to the full theory than the previous one, which we 
illustrate with the following examples in the context of this paper. 

First, the volume operator is constructed immediately from flux operators and does 
not contain functions of connection components which need to be rewritten before they 
can be represented on the Hilbert space. Thus, the volume operator is less ambiguous 
than in the pure connection variables of a feature shared by both the full volume 
operator jSSl IMl 111] and that of homogeneous models jHII. In this context we can also 
discuss the issue of level splitting which was observed even for a single vertex contribution 
of the volume operator in [3T]. With the simple volume operator derived here such a level 
splitting does not occur, and the vertex spectrum is identical to that of a homogeneous 

^This is certainly complicated in loop quantum gravity, not the least because the specific full theory 
is unknown so far. Nevertheless, if models are widely studied and general realizations with welcome 
properties, such as those concerning singularities or conceptual and phenomcnological aspects, have been 
identified, one should hope that a full theory can be realized in such a manner that it reduces to those model 
situations in corresponding regimes. This viewpoint puts the burden on constructing the full theory as well 
as understanding its reduction. The issue is general since one is always forced to employ approximations, 
whether by reduction or otherwise, to understand physical applications of the full theory. Without those 
applications, blind constructions of possible full theories are physically empty. 
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model which has a single rotational axis through each point. Level splitting then occurs 
only when several vertices are considered, which explicitly brings in the inhomogeneous 
properties. 

Secondly, the construction of the Hamiltonian constraint operator fits in the general 
scheme developed in the full theory in [HH] and generalized to homogeneous models in 
|56llHllU]. In particular the presence of possibly non- vanishing covariant components of the 
spin connection in a symmetric model has to be dealt with in a special way compared to the 
full theory. This at first makes it more difficult to interpret the results since the relation to 
the full theory is not as close. However, the fact that it is possible to treat all homogeneous 
models in the same way strongly suggests that there is a general procedure and results are 
not caused by artefacts in a symmetric context. In this paper we have seen that this general 
procedure even extends to inhomogeneous symmetric models which further supports the 
whole construction. The new variables introduced here, even though they were motivated 
independently by implying a simplified volume operator, automatically implement the 
required subtraction of the spin connection in homogeneous components. 

This is a very non-trivial test of the loop quantization procedure: by using these vari- 
ables both the volume operator and the Hamiltonian constraint become closer to what is 
known from homogeneous models and the full theory. Moreover, this works only when the 
special form of spherically symmetric spin connections and extrinsic curvatures (or those 
of the Einstein-Rosen model for which the same procedure works) is taken into account. 
Otherwise, it would be far from clear that the canonical transformation employed here to 
make triad components into momenta amounts to a subtraction of homogeneous spin con- 
nection components. The spherically symmetric model is in between homogeneous models 
and the full theory, and indeed one can observe both homogeneous and inhomogeneous 
aspects. The subtraction is done only in the angular components which correspond to 
homogeneous directions, while for inhomogeneous directions we still have to use the con- 
nection component Ax as would be the case in the full theory. This automatically arises 
from the canonical transformation, and is essential in obtaining a constraint operator with 
the correct classical limit. 

One may think that by way of our canonical transformation we go back from Ashtekar 
variables to ADM like variables with extrinsic curvature as configuration variable. This is, 
however, only partially true since we do not use the x-component of extrinsic curvature, 
but the Ashtekar connection component instead. Thus, only homogeneous directions are 
affected by the transformation, while the inhomogenous components are left unchanged. In 
fact, using extrinsic curvature components throughout, which would be possible from the 
point of view of the symplectic structure, would result in a different Hamiltonian constraint 
operator. For instance, the dependence on t] in (jl^ results from insertions of A in 
not from rj in holonomies as it would happen if Ax was replaced by Kx- Retaining Ax as 
configuration variable and as an argument of basic holonomies is crucial for an operator 
constructed along the lines of the full theory. It is certainly possible to construct constraint 
operators also with ADM like variables, possibly after employing a Bohr representation as 
suggested recently But such a quantization, in contrast to one following a general 

scheme for the full theory and symmetric models as employed here, does not have contact 
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with full quantum gravity which makes reliable conclusions more problematic. 

Moreover, in ADM variables exponentials are unmotivated in contrast to using holonomies 
in theories of connections. As demonstrated by different examples in l53| IT2] , a repre- 
sentation with discontinuous exponentiated connection or extrinsic curvature components 
is induced in models of loop quantum gravity by the full representation. In jSI] it is also 
shown that in one-dimensional models one could construct many different representations 
where also fluxes could be quantized by discontinuous exponentials while the conjugate 
connection components have direct actions without exponentiation. The representation is 
thus far from unique if only the model is considered, and physical properties can crucially 
depend on it. Only by relating models to the full theory with its unique representation 
[nni HOI can a reliable basis for using one representation be given. From this perspective, it 
may be worthwhile to view the representation constructed in [71j, which is related to that 
of jSHl by a gauge fixing, also as a gauge fixing done in a loop formulation removing the 
connection components which from our point of view cannot be replaced by extrinsic 
curvature components. This could provide an embedding of that representation in loop 
quantum gravity. 

6.2 Key simplifications in symmetric models 

While the general scheme discussed in the preceding subsection shows that results can 
be trusted as general loop properties, rather than properties belonging to a given model 
which may or may not be shared by other models or a full theory, it also allows one to see 
how key simplifications arise. They have already been exploited in physical applications 
of homogeneous models, and are now also realized in diagonal inhomogeneous models 
such as the spherically symmetric one or Einstein-Rosen waves. With these models a 
large class of physical situations, which already has been intensively studied classically 
and with diverse quantum methods, is now accessible to explicit computations in the loop 
framework. Applications span all areas of gravitational physics from cosmology to black 
holes and gravitational waves. 

A technical feature shared by all these models is that the relevant form of invariant 
connections and triads is diagonal, i.e. Lie algebra valued components corresponding to 
independent directions such as x, ip here are perpendicular in the internal direction. 
How this leads to simplifications has been explained in detail in ^Tj. This diagonalization 
essentially amounts to a reduction from SU(2) to U(l) which is the reason for a simpler 
volume operator. However, as seen by comparing this paper with (321; this reduction in 
the gauge group is necessary but not sufficient for simplifying the calculations. Also an 
appropriate form of canonical variables is essential, which connects simplifications in the 
volume operator with the construction of the Hamiltonian constraint. 

The required computations then are quite similar in homogeneous and inhomogeneous 
models. Only the discussion and solution of the constraint equation is, of course, more 
complicated in inhomogeneous situations since more degrees of freedom are involved. But 
in all these models the constraint operator takes an analogous form, and in all models 
a triad representation is available. The latter feature is the main difference to the full 
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theory, besides the comphcation in exphcit calculations. Writing the constraint equation 
in the triad representation is much more intuitive than in the connection representation 
and has led to many results concerning the quantum structure of classical singularities 
[milHlIZaiZSllZai^linilSliniEn!- Again, one may worry that this difference to the full 
theory implies special features of models which may be misleading. However, this is not 
the case since we just transform an equation which we have in the models as well as in 
the full theory into a different form more suitable to applications. We are, however, not 
changing the equation itself or its solutions. 

In this triad representation we then obtain difference equations, which are of different 
complexity: ordinary in the vacuum isotropic case, partial in homogeneous models or even 
partial in infinitely many variables in inhomogeneous models. In all cases, the nature of 
being a difference equation is derived from holonomy operators appearing in the constraint, 
with coefficients determined by commutators with the volume operator. Computing those 
coefficients can be done with equal ease and no new ingredients in all models, and is crucial 
for the discussion of singularities (see also for a recent summary). 

Another shared feature of the models is the relation between the Euclidean and Lorentzian 
constraints, which is always simple and only involves a rescaling of the holonomy part of 
the constraint. This is in contrast to the full theory where the relation is much more 
complicated and the Lorentzian constraint, so far, can only be quantized by introducing 
additional commutators between the Euclidean constraint and the volume operator [HH] . 
At this point it is not clear if the simplification in symmetric models is general enough to 
lead to reliable results. This can, however, be tested directly since the procedure of the full 
theory is certainly available in the models, too. Usually, one first studies the simpler ver- 
sion in order to understand the physical implications, before more complicated choices can 
be compared (see, e.g., [71] compared to jZH]). So far, the more complicated version only 
implies higher order equations and more involved coefficients, but no crucial differences. 

6.3 Physical applications 

The key simplifications discussed before indicate that physical applications can be obtained 
in an explicit form even at the dynamical level. There are several features common with 
the formalism of homogeneous models, such as the explicitly known matrix elements of 
the Hamiltonian constraint in the Euclidean and Lorentzian form, the fact that they are 
sufficiently simple functions of the spin network labels and, crucially, the availability of 
a triad representation. The latter facilitates in particular the analysis of a neighborhood 
of classical singularities in order to see if they still present a boundary to the quantum 
evolution. This is crucial for an interpretation, not the formulation of models; the lack of 
a triad representation in the full theory is thus not problematic for trusting the model. 

However, there are certainly also complications as compared to homogeneous mod- 
els with only finitely many kinematical degrees of freedom. The spherically symmetric 
quantum constraint equation in the triad representation is a partial difference equation in 
potentially infinitely many variables and, depending on the form of the quantization, the 
number of degrees of freedom involved at each time step may not even be constant (if the 
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constraint operator creates new vertices). Moreover, there are many more ways to approach 
a classical singularity on midisuperspace, depending also on gauge choices, and there are 
different versions of singularities. This makes general statements about their removal or 
persistence more complicated since for that classical singularities first have to be located 
in midisuperspace before the possibility of unique extensions of physical wave functions 
is discussed. Yet, in one-dimensional models this is possible and shows many crucial new 
features The ordering, for instance is more restricted than in homogeneous models 
and a symmetric one is preferred. While homogeneous models are non-singular for the 
non-symmetric and symmetric ordering because the evolution can be continued even in 
cases where coefficients of the difference equation can become zero ^\ , in inhomogeneous 
models this is not possible at least for the direct operator following from a non-symmetric 
version. As noted before, the symmetric ordering, on the other hand, has generically non- 
zero coefficients in such a manner that the evolution does not break down. Requiring a 
non-singular evolution thus restricts the ordering choices in inhomogeneous models. 

For other physical applications it is helpful to have approximation schemes available 
which allow one to isolate the essential quantum modifications to classical equations. In 
homogeneous models, effective classical equations of the form introduced in [121 113 HH] 
have been essential in many recent physical applications in the context of cosmology fK]\ 
nii^M^nBE^^E^^^^E^EM^EB- These are classical equations in that 
they are ordinary differential equations in coordinate time which are much easier to handle 
than difference equations [IHll^lZlllZSllSnilHZllHHllHniinnilin!. Quantum effects are then 
imported by comparing the expectation value of the quantum constraint with the classical 
constraint equation, which can also be analyzed and justified by comparing the motion of 
quantum wave packets with solutions to the effective equations OIlSj. (This is part of a 
general procedure which also contains the usual effective action techniques jHH 1221 ESI ES]-) 
Similarly, effective classical equations for inhomogeneous models, which then would be 
partial differential equations in space-time coordinates, would be much easier to deal with 
than a difference equation in infinitely many variables. Moreover, mathematical expertise 
from geometrical analysis would be available to arrive at general conclusions. 

In cosmological models, several quantum effects have been observed which give rise 
to different terms in effective classical equations. The first and most direct one was the 
modification of matter Hamiltonians at small volume |S3J ESI, which makes classically 
diverging energy densities finite along effective trajectories and also played a role for the 
removal of cosmological singularities. The same modification is available in a spherically 
symmetric model with matter; but since here even a vacuum model would have a classical 
(Schwarzschild) singularity, it cannot be expected to be sufficient for a non-singular effective 
formulation. 

Other modifications come from the gravitational part of the constraint and are thus 
available even in the vacuum case. There are perturbative corrections which can be in- 
terpreted as being analogous to higher curvature terms, and non-perturbative ones in the 
case of a non-vanishing spin connection. (The latter can be viewed as providing a natural 
cut-off on intrinsic curvature, analogously to the extrinsic curvature cut-off from a quan- 
tum matter Hamiltonian in a homogeneous situation.) The non-perturbative modifications 
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have mainly been used so far in the Bianchi IX model where they remove the classical chaos 
[221IH3- They are derived by quantizing the spin connection potential term of the Hamil- 
tonian constraint in accordance with the general scheme described above. Since the spin 
connection contains inverse powers of the triad components, it will obtain modifications 
from the quantization of inverse powers IHK] similarly to a matter Hamiltonian. 

Here, the same procedure is available where we also have an intrinsic curvature po- 
tential with the spin connection depending on inverse powers of E'^ . Using eigenvalues of 
the quantized potential in the classical constraint equation then provides effective classical 
equations showing one main quantum effect from the cut-off of intrinsic curvature. Com- 
parison with the Schwarzschild solution shows that this indeed provides modifications at 
the right place, namely close to the classical singularity where E"^ is small. In asymptotic 
regimes or around the horizon of massive black holes, on the other hand, E"^ is large and 
so the behavior there remains classical. This would be different had we used other metric 
variables such as the co-triad or metric instead of the densitized triad as, e.g., in pH] . 
Then the component whose inverse appears in the spin connection (|2U|) would be e^-, which 
for Schwarzschild is large at the horizon but also at the classical singularity and so there 
would be no quantum corrections there. In asymptotic regimes, on the other hand, 
would be of the order one: modifications would be noticable but completely unwanted. 
Thus, similarly to homogeneous models jH] we see that the issue of singularity removal will 
crucially depend on the canonical variables used. Densitized triad components, which we 
have to use anyway since they are part of the basis of the full background independent 
quantization, are well-suited for this aim. These statements are only preliminary as of now 
because no consistent set of effective equations for inhomogeneous cases has been derived 
so far. 

Besides singularities, horizons would be the second interesting aspect to be studied. 
As just noted, we do not expect strong quantum modifications there at least for massive 
black holes. But there are other interesting aspects which are usually related to quantum 
models, such as the issue of black hole entropy, horizon degrees of freedom and fluctuations, 
and Hawking radiation. All these issues, for instance in the context of the scenario of 
|27j . require solving the Hamiltonian constraint which can provide important feedback on 
the viability of a given quantization scheme, and thus test the full theory and restrict 
quantization ambiguities. The isolated or dynamical horizon framework jHH IHS HEl EZl 
EHl inni provides an ideal setting for an analysis at the classical and quantum levels. If 
a horizon is isolated or slowly evolving jlUUj . certain terms in the quantum constraint 
become negligible such that the constraint simplifies |lUlj . This provides an approximation 
scheme to understand the physics at the horizon or, in a perturbative form, around the 
horizon. Also dynamical processes can then be studied in a controlled manner because 
not just isolated but also slowly evolving horizons are allowed. In this way one can derive 
physical information about black hole or other horizons, but also study general issues of 
the Hamiltonian constraint such as observables. 
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